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ABSTRACT 

We demonstrate the ability to measure precise stellar barycentric radial velocities with the dispersed 
fixed-delay interferometer technique using the Exoplanet Tracker (ET), an instrument primarily designed 
for precision differential Dopplcr velocity measurements using this technique. Our barycentric radial 
velocities, derived from observations taken at the KPNO 2.1 meter telescope, differ from those of Nidever 
et al. by 0.047 km s -1 (rms) when simultaneous iodine calibration is used, and by 0.120 km s" 1 (rms) 
without simultaneous iodine calibration. Our results effectively show that a Michelson interferometer 
coupled to a spectrograph allows precise measurements of barycentric radial velocities even at a modest 
spectral resolution of R ~ 5100. A multi-object version of the ET instrument capable of observing 
~500 stars per night is being used at the Sloan 2.5 m telescope at Apache Point Observatory for the 
Multi-object APO Radial Velocity Exoplanet Large-area Survey (MARVELS), a wide-field radial velocity 
survey for extrasolar planets around TYCHO-2 stars in the magnitude range 7.6 < V < 12. In addition 
to precise differential velocities, this survey will also yield precise barycentric radial velocities for many 
thousands of stars using the data analysis techniques reported here. Such a large kinematic survey 
at high velocity precision will be useful in identifying the signature of accretion events in the Milky 
Way and understanding local stellar kinematics in addition to discovering exoplanets, brown dwarfs and 
spectroscopic binaries. 

Subject headings: techniques: radial velocities — techniques: spectroscopic — instrumentation: 

interferometers — instrumentation: spectrographs — stars: kinematics — methods: 
data analysis 



1. INTRODUCTION 

The discovery of various tidal streams in our galaxy 
(Ibata et al. 1994; Yanny et al. 2006; Belokurov et al. 
2006) challenges the monolithic collapse formation model 
(Eggen 1962), but the streams discovered so far form only 
a small fraction of the mass of the galaxy. Large kinematic 
surveys are necessary to determine if such merger and ac- 
cretion events are more common. Radial velocity surveys 
using multi-object spectrographs, like RAVE (Steinmetz et 
al. 2006), are already producing barycentric radial veloci- 
ties for many thousands of stars at a precision of 2 — 3 km 
s _1 , enabling studies of galactic dynamics (e.g. Veltz et 
al. 2008). In the northern hemisphere the magnitude lim- 
ited Geneva-Copenhagen Survey (Nordstrom et al. 2004) 
has provided radial velocities for over 14000 bright nearby 
F and G dwarfs (V < 8.6) at a precision of 100 — 300 m 
s _1 . Coupled with gravity and chemical composition de- 
rived from the observed spectra, such datasets of precise 
barycentric radial velocities are very valuable in identify- 
ing past accretion events and substructure in the galac- 
tic disk (Helmi et al. 2006), finding members of moving 
groups (Jones 1970) and determining cluster membership. 

The ability to determine precise differential radial veloc- 
ities has also improved dramatically in the last few years, 
and planet search surveys are currently at a precision of 
1 - 3 m s" 1 (Rupprecht et al. 2004; Butler et al. 1996). 
Such planet search surveys, however, use an arbitrary zero 



point for the stellar velocities since only differential veloc- 
ities are measured. The true radial velocity of the center- 
of-mass of a star is a difficult quantity to measure accu- 
rately using spectroscopy since it involves additional mod- 
elling (or assumptions) of effects like convective blue shift 
and gravitational redshift (Lindegren & Dravins 2003). A 
quantity that is easier to obtain, in terms of observables 
and known physical constants, is the apparent Doppler 
shift of the lines in the stellar spectra, as measured in the 
solar system barycenter, with the Sun setting the veloc- 
ity zero point. This Doppler shift can be converted into a 
velocity that we define as the barycentric radial velocity 
of the star. Very precise barycentric radial velocities have 
been measured using high resolution echelle spectrographs 
by Udry et al. 1999 (UD99 hereafter) for a number of 
stable F, G, K dwarfs, as well as by Nidever et al. 2002 
(ND02 hereafter) for a larger sample. 

In this paper we demonstrate the ability of a dis- 
persed fixed-delay interferometer (DFDI) to measure pre- 
cise barycentric radial velocities. The Exoplanet Tracker 
(ET) is a prototype of this class of instruments, which 
use a Michelson interferometer followed by a medium res- 
olution spectrograph working in the first grating order. 
Fixed delay interferometers have long been used by so- 
lar astrophysicists to measure solar oscillations (Gorskii 
& Lebedev 1977; Kozhevatov 1983; Harvey et al. 1995). 
In 1997 D. J. Erskine proposed adding a spectrograph as 
a post disperser for increasing the measurement precision 
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of differential radial velocities (Erskine & Ge 2000; Ge, 
Erskine & Rushford 2000). When using this technique si- 
nusoidal interference fringes are formed in the slit direction 
at the position of the stellar absorption lines. A shift in 
the dispersion direction of the underlying spectra due to 
a Dopplcr shift manifests itself as the phase shift of the 
sinusoidal fringes in the slit direction (Ge 2002). 

Precise radial velocities can be measured with DFDI 
even at low wavelength dispersions since information is 
also being encoded as sinusoidal fringes in the slit direc- 
tion. This results in each stellar spectrum taking up only 
a small fraction of the available number of pixels on the 
CCD detector used to record the spectrum. This feature 
of the DFDI technique is a great advantage since it has 
the potential to enable a multi-object instrument (similar 
to ET) to simultaneously obtain precise radial velocities 
for a large number of stars (Ge 2002; Ge et al. 2002) with- 
out the need for an expensive cross-dispersed high spectral 
resolution instrument. One such instrument, the W. M. 
Keck Exoplanet Tracker, is being used at the Sloan 2.5 m 
telescope for the magnitude limited Multi-object APO Ra- 
dial Velocity Exoplanet Large-area Survey (MARVELS), 
a survey for extrasolar planets around TYCHO-2 stars in 
the magnitude range 7.6 < V < 12 that is being con- 
ducted as part of SDSS III (Weinberg et al. 2007). This 
instrument can observe 59 stars simultaneously and ap- 
proximately 500 stars per night. Designed primarily as a 
wide-field precision radial velocity survey to discover short 
and intermediate period planets, MARVELS will also pro- 
vide a dataset of precise barycentric radial velocities for 
thousands of late F, G, and K TYCHO-2 stars. 

2. OBSERVATIONS 

The observations presented in this article were con- 
ducted with the single object ET instrument at the Kitt 
Peak 2.1 meter telescope in January 2006 as part of an 
ongoing planet search program. The ET instrument is 
linked to the telescope by a 200 /xm optical fiber, which is 
fed by the f/8 beam of the 2.1 meter telescope. The instru- 
ment itself is housed in an isolated room in the basement 
and consists of a Michelson interferometer in series with 
a medium resolution spectrograph operating in the wave- 
length range 5000 — 5640 A, and with a spectral resolution 
of R ~ 5100. An iodine cell can be inserted into the optical 
path and serves to calibrate out instrument drifts. This 
instrument has successfully demonstrated its short-term 
stability by confirming known planets (e.g. van Eyken et 
al. 2004) and discovering a short period planet around the 
star HD102195 (Ge et al. 2006). 

The observing procedure for the ET planet search pro- 
gram is similar to other precision radial velocity programs 
that use an iodine cell to calibrate out the instrument drifts 
(e.g. Butler et al. 1996). For each star in the program a 
single exposure is acquired without the iodine cell inserted 
in the beam path. An exposure of the iodine cell, illumi- 
nated by a tungsten lamp, is also acquired immediately 
before or after this stellar exposure. These exposures are 
referred to as the star and iodine template respectively. 
The iodine template is acquired at high signal to noise 
(S/N), typically S/N = 150-200 per pixel. The signal-to- 
noise on the star template depends on the brightness of the 
star and the exposure time and is generally in the range 



of S/N = 50-200 per pixel. All subsequent exposures are 
taken with the iodine cell inserted in the path of the stel- 
lar beam. It is from these exposures that precision radial 
velocities are obtained since the velocity drift of the iodine 
component tracks the instrument drift, allowing it be cal- 
ibrated out. Exposure times are similar to those used on 
the stellar template, typically yielding a S/N ~ 1/V2 of 
that obtained on the stellar template (since the iodine cell 
absorbs almost 50% of the light). 

Stable stars and known planet-bearing stars are rou- 
tinely observed for calibration as part of the ET planet 
search program. During the observing run we were able 
to obtain regular observations of the RV stable stars -q 
Cas (small linear trend in velocity), r Ceti, and 36 Uma 
as well as the known planet bearing stars 51 Peg (Mayor 
& Queloz 1995), 55 Cnc (McArthur et al. 2004), and 
HD4967 (Butler et al. 2003). For six stars, chosen be- 
cause they had rms velocity scatter less than 100 m s^ 1 re- 
ported in ND02 (HD3674, HD3861, HD12414, HD105405, 
HD130087, HD134044), we also acquired stellar templates 
and iodine templates immediately before or after the re- 
spective stellar templates, although no star-with-iodine 
data points were obtained for these six stars. Typical ex- 
posure times varied from 3 minutes for the brightest star (77 
Cas, V = 3.45) to 10 minutes for the faintest (HD130087, 
V = 7.52), leading to a photon noise limited radial velocity 
error of 3 — 10 m s _1 for most observations. 

The advantage of an isolated, fiber-fed instrument is the 
potential for high stability. The total instrument velocity 
drift during the observing run (as tracked by the iodine 
templates) can be seen in Figure 1. The instrument ve- 
locity drift is only 2.5 km s _1 over the entire 12 day ob- 
serving run and less than 0.5 km s _1 over the first night. 
This short term stability is essential for the data analysis 
technique described in Section 3.1. 

3. DATA ANALYSIS PROCEDURE 

The standard ET data analysis pipeline was designed 
to extract differential radial velocity for the planet search 
program. It does not give the barycentric radial velocity of 
the star since the zero point itself is arbitrary. The pipeline 
processing steps are described in more detail in van Eyken 
et al. 2004, and Ge et al. 2006. Here we discuss the 
various processing steps briefly, but we concentrate on the 
modifications that allow barycentric radial velocities to be 
measured. 

The data produced by ET were processed using stan- 
dard IRAF procedures, as well as software written in Re- 
search System Inc.'s IDL software. The images were cor- 
rected for biases, dark current, and scattered light and 
then trimmed, illumination corrected, and low-pass fil- 
tered. The visibilities (V) and the phases (9) of the fringes 
were determined for each channel by fitting a sine wave to 
each column of pixels in the slit direction, with the fringe 
visibility, V, being defined as 

V {I'max 1 rain) / '(^ max ~t~ Imin) (1) 

where I ma x and I m in are the maximum and minimum val- 
ues of the fitted sine wave. A wavelength calibration was 
applied using the observed Thorium- Argon spectra, and 
the visibility and phase data were dispersion corrected 
and re-binned to a log-lambda wavelength scale. To de- 
termine differential velocity shifts in the star+iodine data 
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(starlight passing through the iodine cell) the data can, 
to a good approximation, be considered as a summation 
of the complex visibilities (V = Ve %e ) of the relevant star 
(Vse l6s o) and iodine (Vie l6 'o) templates (van Eyken et al. 
2004; Erskinc 2003). For small velocity shifts the complex 
visibility of the data can be written as 

V D e i0D = V s e ies oe i0s - ies ° + Vie™ 1 * e ie '~ ie '° (2) 
where Vd,Vs, and Vj are the fringe visibilities for a given 
wavelength in the star+iodine data, star template and io- 
dine template respectively, and 8d,6s , and 9j the corre- 
sponding measured phases. In the presence of real velocity 
shift of the star and instrument drifts the complex visibil- 
ities of the star and iodine template best match the data 
with a phase shift of 9s — 9s and 9j — 6j respectively. 
The iodine is a stable reference and the iodine phase shift 
tracks the instrument drift. The difference between star 
and iodine shifts is the real phase shift of the star, A<f>, 
corrected for any instrumental drifts 

AJ> = (9 S - 9 So ) - (0j - 0/ o ) (3) 
This phase shift can be converted to a velocity shift (Av) 
by a known phase-to- velocity scaling factor that is a func- 
tion of the optical delay in the Michelson interferometer 
(d), the wavelength (A) and the speed of light (c). These 
parameters are connected by the relation 

In reality the phase shifts of the star and iodine tem- 
plates are not the only free parameters needed to fit the 
observed data. The Doppler velocity shift of the star in- 
duces a wavelength shift of the entire spectrum and in- 
strument instabilities can lead to pixel shifts relative to 
the templates. Finding the correct phase shifts requires 
that these wavelength shifts (or pixels shifts) for the star 
and iodine templates also be fit simultaneously with the 
phase shift parameters. 

To determine barycentric radial velocities we used the 
stellar template of a star with known barycentric radial ve- 
locity (n Cas) in place of the target star template in Equa- 
tion 2. The visibility of the reference star was allowed to 
scale to partially compensate for the template mismatch. 
The stellar visibility depends on the line depth and width 
(Ge 2002), and allowing it to scale is in effect similar to 
attempting to match the line depths and widths of the 
stellar template with that of the observed star. Similar 
spectral-morphing techniques have been used (Johnson et 
al. 2006) with high-resolution echelle data from the N2K 
survey to achieve short term differential velocity precision 
of ~ 5 m s" 1 without the need for a stellar template. 

3.1. Radial Velocity using Templates Alone 

The stability of the ET instrument (Figure 1) makes it 
possible to attempt to extract the radial velocity without 
using simultaneous calibration with the iodine cell. When 
the cell is not inserted into the optical path, the resulting 
observations of the star can explicitly be modelled as the 
template spectrum with a phase shift added (the template 
spectrum may also be the stellar template of a reference 
star) 

V D e l0D =v s e if)s °e i9s - ies ° (5) 

In such a case the phase shift, 9s — 9s , is a combination 
of the velocity shift of the star compared to the template 



and the instrument drift. When using this technique the 
ET observing procedure requires all star-only exposure to 
have an iodine template taken immediately before or after 
the exposure. If the instrument is stable enough, these io- 
dine exposures can be used to determine and calibrate out 
the instrument drifts for the star exposures. This is only 
possible when the instrumental velocity drifts are small in 
the time interval between the exposure mid-points of the 
star and iodine templates. 

Using the entire spectral range of ~ 640 A a cross- 
correlation was performed in log-lambda space using the 
visibilities (V) of the reference star and the target star. 
The purpose of this cross-correlation was to yield an ap- 
proximate velocity shift, accurate to ~ 1 — 2 km s -1 , which 
is useful as a first guess starting point to calculate a more 
precise velocity shift using both the visibility and phase 
information. 

The velocity extraction algorithm that is part of the ET 
pipeline was used to minimize the residuals between the 
complex visibilities (V) of the reference star and the tar- 
get star by varying the phase shift applied to the reference 
template complex visibility. The visibility of the reference 
star was also allowed to scale in the fitting procedure in an 
attempt to minimize the template mismatch, as described 
in the previous section. This technique yields a phase shift 
&s — ds that is a measure of the difference in velocity be- 
tween the stellar template used and the target star. 

To account for instrument drifts between the reference 
star template and the target star, a similar analysis pro- 
cedure was followed for the iodine templates (taken very 
close in time to the reference and target star exposure) 
with the only difference being that the iodine visibility was 
not allowed to scale (since there is no template mismatch 
between iodine cell spectra taken at different times). The 
phase drift between the two iodine templates is a direct 
measure of the instrument drift. Subtracting off the io- 
dine phase shift corrects for the instrument drift, giving 
the real phase shift of the target star with respect to the 
reference star. 

3.2. Radial Velocity from Star+iodine Data 

In the normal data acquisition mode the stellar beam 
passes through the iodine cell. Cross-correlation with the 
visibilities from a reference stellar template alone is not 
appropriate since the iodine spectrum is also present in 
the data. Using Equation 2 we can write 

Vd = Vl + V? + 2VsVj cos(9 s - 0j) (6) 
The term 2VsVi cos(#s — #r) has an expectation value 
close to zero over the 640 A wavelength band since the 
cosine term fluctuates randomly. We ignored this term 
and approximated the visibilities as adding in quadrature 
(Vg ~ y§ + Vj 2 ). This allowed us to use the TODCOR 
algorithm (Zucker & Mazeh 1994) on the ensemble of mea- 
sured visibilities to determine the approximate velocity 
shifts of the star and iodine templates that are a best fit to 
the data. The TODCOR algorithm is a two dimensional 
cross-correlation technique generally used to determine the 
individual radial velocities of both components of a spec- 
troscopic binary, and can be applied in this case since the 
two components are approximately additive. Since a refer- 
ence stellar template for a different star (n Cas) was being 
used, a visibility scaling was derived using mean visibilities 
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and Equation 6, and used to scale the stellar visibilities for 
TODCOR. The scaling factor in visibilities accounts par- 
tially for the mismatch in stellar templates. The approxi- 
mate radial velocity values derived using TODCOR were 
useful in determining the correct phase shift when the shift 
was more than ±2ir. These radial velocities were used as 
a first guess to determine the correct wavelength-shift and 
phase-shift of the stellar and iodine templates that mini- 
mized the residuals in the complex visibility between the 
data and the best fit solution. Subtracting off the iodine 
phase shift gives the phase shift of the target star with re- 
spect to the reference star, corrected for instrument drift. 

4. RADIAL VELOCITY RESULTS 

The differential radial velocity precision of one of the ref- 
erence stars, r\ Cas, is shown in Figure 2. These velocities 
were extracted using the techniques described in Section 
3, but with an arbitrary zero-point. In the case of r\ Cas 
the velocity results are identical to those obtained using 
the standard pipeline since this object is the reference star. 
The velocity rms of ~ 6.4 m s _1 is due to a combination of 
photon noise errors in the data and the templates and sys- 
tematic errors in the velocity extraction. We used this star 
as the stellar reference template to determine barycentric 
radial velocities. Velocities for all the stars are calculated 
relative to n Cas, requiring prior knowledge the barycentric 
radial velocity of n Cas itself to apply the correct velocity 
offset. We used the r\ Cas mean barycentric radial velocity 
from ND02 (8.314 km s _1 ). 77 Cas is known to be a long- 
period binary, whose velocity is well fit with a linear trend. 
We transformed the ND02 velocity to our epoch of obser- 
vation by applying a correction term for its known 9.3 m 
s _1 year -1 linear velocity drift (Cummings ct al. 1999). 
Any error in the determination of this velocity is a error 
in the zero-point itself, and will affect all target stars in 
a similar way. The phase difference of each star from the 
77 Cas stellar template was converted to velocity using the 
phase-to- velocity scaling parameter (c\/2nd). This radial 
velocity was then corrected to the solar system barycenter 
by removing the barycentric velocity contribution of the 
telescope, calculated with software written in IDL that 
use the JPL ephemeris. These algorithms are accurate to 
better than ~ 2 m s _1 . The corrected barycentric radial 
velocity of 77 Cas was added to obtain the barycentric ra- 
dial velocity of the star. 

4.1. Radial Velocity Results from Templates Alone 

Barycentric radial velocities were derived using the pro- 
cedure outlined in Section 3.1. A comparison of these ve- 
locities with those published by ND02 is shown in Figure 
3. No attempt has been made to correct the velocities 
for the known planetary companions to 51 Peg, 55Cnc 
and HD49674, which may cause additional velocity scat- 
ter. Our velocities agree well with those of ND02 and show 
a scatter of 120 m s _1 , with a mean offset of —79 ± 25 m 
s _1 . The agreement of our results with those of ND02 con- 
firms the validity of using the iodine templates to remove 
the instrument drift when such drifts are small (Figure 
1). For example, the typical instrument drift over the first 
night is approximately 10 m s _1 over a 15-minute expo- 
sure. Such exposure times with the ET instrument allow 
the radial velocity of an II th magnitude (V=ll) star to 



be determined at this level of precision (100 m s 1 ). 

4.2. Radial Velocity Results from Star+Iodine Data 

When the required precision is high, it is perhaps more 
appropriate to use the star+iodine data since the iodine 
tracks the instrument drift exactly. For all the reference 
and planet bearing stars, the barycentric radial velocities 
obtained for all data points (using 77 Cas as template), 
derived using the procedure outlined in Section 3.2, were 
averaged together to give a mean barycentric radial ve- 
locity for each star. These radial velocities are shown in 
Figure 4 and display a rms scatter of 47 m s _1 with a mean 
offset of -3 ± 19 m s _1 in comparison with ND02. 

4.3. Systematic Errors and Velocity Zero Point 

Tests with twilight solar spectra as the reference stellar 
template show similar levels of rms scatter to those using 
77 Cas, but the twilight observation of the solar spectrum 
(instead of day sky) introduced an offset in the velocities 
when the solar spectrum was used as the template. Our 
results depend on the velocity of 77 Cas derived by ND02 
who use the National Solar Observatory (NSO) solar spec- 
trum (Kurucz et al. 1984) as the stellar template with a 
522 m s _1 offset to make the radial velocity of the Sun 
zero. Due to the use of the same zero point our estimated 
velocities will exhibit similar systematic errors with spec- 
tral type as ND02, in addition to any template mismatch 
effects introduced due to our use of n Cas (GOV) instead 
of the Sun (G2V). 

A number of systematic errors can occur in velocity de- 
termination using this technique. Fast rotators (wsinz > 
12 km s _1 ) have very low fringe visibilities at our chosen 
interferometer optical delay (7 mm), leading to a signifi- 
cantly reduced velocity precision for such stars. Incorrect 
calibration of the interferometer delay can lead to errors 
in the phase-to-velocity scaling, causing errors in the cal- 
culated radial velocity. Our reference (77 Cas) has a small 
linear drift in velocity that we have accounted for using 
the mean epoch of observation in ND02. An error in this 
determined barycentric radial velocity of 77 Cas will appear 
as a shift in the mean offsets reported here. The stars with 
known planets have additional velocity variability due to 
their companions. This variability is reduced when using 
star+iodine data (Section 4.2) and averaging all the ob- 
served velocities, especially for the short period planets 51 
Peg and HD49674 that have orbital periods less than our 
12 day observing run. The multi-planet 55 Cnc system 
(McArthur et al. 2004), though, has long period planets 
which affect the calculated radial velocity. For each star in 
their sample ND02 report only the mean of the observed 
barycentric radial velocities, and to ensure a consistent 
comparison we have made no attempt here to correct our 
data for these known planetary companions. In general, 
velocity corrections for known companions is easy to im- 
plement if the orbital parameters are known, but we have 
not done so for consistency with ND02. 

5. DISCUSSION 

We have demonstrated the ability to measure barycen- 
tric radial velocities at a high level of precision (rms = 
50 — 120 m s _1 ) using a dispersed fixed delay interfer- 
ometer instrument. Our results also show that the radial 
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velocities derived using a very different technique and an 
instrument at substantially lower resolution are not very 
different in precision from those derived using high resolu- 
tion spectroscopy and cross-correlation with stellar (ND02; 
UD99) or Fe I line templates (Gullberg & Lindegren 2002). 
Simultaneous calibration with the iodine cell yields more 
precise velocities, but the data analysis using only the tem- 
plates (Section 3.1, 4.1) may be more appropriate for faint 
stars where the low signal to noise may introduce signifi- 
cant errors when extracting velocities from the star+iodine 
data. We have corrected the observed radial velocities to 
the barycenter of the solar system. Using the corrected rj 
Cas radial velocity from ND02 to set the velocity scale for 
our radial velocities then automatically corrects all veloc- 
ities for the gravitational redshift and convective blueshift 
of the Sun (used by ND02 as the zero point) . By adopting 
the radial velocity of rj Cas from ND02 we have automati- 
cally adopted ND02 as our velocity reference for the abso- 
lute calibration of a zero-point and performed a secondary 
measurement to determine the velocities of all our other 
target stars. No attempt is made to further correct the 
stellar radial velocities for spectral type and luminosity 
class dependent effects. This results in our derived veloci- 
ties having typical accuracies of 0.2 — 0.3 km s _1 (derived 
from ND02, our velocity reference) for late F, G and K 
dwarfs, and precisions of ~ 50 — 100 m s _1 when compar- 
ing stars of similar spectral types. With the ET instrument 
scheduled to be a facility instrument at the KPNO 2.1m, 
this velocity measurement capability will be available to 
the wider astronomical community. These techniques will 



also be applied to data from the MARVELS survey that 
will observe ~ 10,000 stars as part of SDSS III, leading 
to precise barycentric radial velocities being measured for 
the majority of the stars in the sample. 

Future work with these instruments will focus on using 
the day sky spectra of the Sun as a stellar template to 
determine the velocity zero point, allowing independent 
comparison with the radial velocities of ND02 and UD99 
without having to adopt either of these works as the veloc- 
ity reference, and a better understanding of any zero point 
differences that may arise due to the much lower spectral 
resolution (Dravins & Norlund 1990) of ET compared to 
high resolution echelle spectrographs. 
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Fig. 1. — Velocity drift of the ET instrument. The solid line is the velocity drift as tracked by the iodine templates (filled circles). The 
dashed line tracks the temperature drift of the Michelson interferometer (asterisks). A loose correlation between the temperature drifts and 
the velocity drifts is evident. Instrument stability is essential for obtaining radial velocities if the simultaneous iodine calibration is not used. 
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Fig. 2. — Differential radial velocity observations of the star -q Cas with ET. The rms velocity scatter of ~ 6.4 m s _1 is caused by a 
combination of photon noise errors and systematic errors. This star is used as the stellar template to determine barycentric radial velocities 
(see Section 4). 
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Fig. 3. — Comparison of radial velocities obtained by ET, from stellar and iodine templates alone, with ND02. The comparison of ET 
velocities with ND02 velocities yields an rms velocity scatter of 120 m s~ 1 . 
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Fig. 4. — Residuals of the average ET radial velocities derived from star+iodine data from the velocities of ND02. Residuals show an rms 
velocity scatter of 47 m s _1 . 



